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ABSTRACT 

This  paper  treats  conformal  mapping  as  it  relates  to  the  generation  of  grids 
to  be  used  for  flow  simulation.  Classical  and  contemporary  mapping  methods  are 
discussed  and  compared  in  some  detail.  Numerous  suggestions  are  included  to 
help  the  reader  avoid  common  pitfalls,  and  the  mapping  methods  believed  most 


promising  are  identified. 


INTRODUCTION 

Conformal  napping  is  a versatile  component  of  the  spectrum  of  grid  generation 
techniques,  yet  its  use  is  sometimes  avoided  by  investigators  who  feel  uneasy 
about  its  implementation  due  to  unfamiliarity.  The  object  of  this  paper  is  to 
put  conformal  mapping  into  perspective,  to  discuss  alternate  implementations, 
and  to  induce  the  reader  to  use  conformal  mapping,  when  appropriate,  as  one 
component  of  a grid  generation  process. 

Although  this  paper  discusses  the  generation  of  grids  using  conformal  mapping 
techniques,  the  grids  themselves  are  not  restricted  to  being  conformal  or 
orthogonal.  Algebraic  techniques  are  considered  for  the  latter  stages  of  a 
mapping  sequence,  and  elliptic  techniques  are  considered  for  "filling-in"  the 
interior  of  a grid  once  the  boundary  correspondences  have  been  obtained  through 
a mapping  sequence.  The  emphasis  is  on  grids  on  a surface,  rather  than  through- 
out a volume,  since  a conformal  mapping  is  basically  a surface-to-surface 
correspondence . 

This  paper  starts  with  a discussion  of  the  basic  differences  between 
algebraic,  orthogonal,  and  conformal  coordinate  systems.  Some  implications  of 
these  differences  are  considered  as  they  relate  to  flow  simulation  and  design 
considerations.  A brief  review  of  complex  variable  notation  and  various  types 
of  conformal  mappings  is  followed  by  a discussion  of  the  problem  of  multiple 
valuedness,  which  can  be  the  most  troublesome  aspect  of  conformal  mapping  from 
a computer  implementation  viewpoint.  Various  mapping  techniques  in  current  use 
are  reviewed  and  techniques  useful  for  creating  new  conformal  mappings  are  out- 
lined. The  generation  of  a grid,  once  the  mapping  is  known,  is  discussed. 
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Suitable  reference  material  is  listed  by  problem  type  and  packaged  computa- 
tional tools  are  described.  The  closing  remarks  highlight  techniques  that  are 
believed  to  be  the  most  promising  in  this  evolving  field. 

COORDINATE  TRANSFORMATIONS 

This  section  illustrates  some  of  the  differences  between  general,  independ- 
ent algebraic,  orthogonal,  and  conformal  coordinate  systems. 

A general  two-dimensional  transformation  between  (x,y)  physical  coordinates 
and  (5,n)  computational  coordinates  (where  (£,n)  are  taken  to  be  orthogonal 
throughout  this  paper)  can  be  written  as 

x - x(£,n) 

y - y(£,n) 

or  in  differential  form  in  terms  of  a Jacobian  matrix  as 


dx 

a b] 

d£ 

s 

| ‘ 

dy. 

c dj 

dn 

where,  from  the  chain  rule,  the  four  independent  parameters  a,  b,  c,  and  d are 
given  by 


d ■ u ? . 

Here  dx  and  dy  are  components  of  a differential  vector  dZ  in  the  physical  plane, 
while  d£  and  dn  are  orthogonal  (perpendicular)  components  of  a differential 
vector  dC  in  a computational  plane.  For  a general  algebraic  transformation,  dx 
and  dy  are  not  orthogonal,  and  their  magnitudes  need  not  be  equal  when  the 
magnitudes  of  d£  and  dn  are  equal.  An  elliptic  technique  as  in  Ref.  [1],  for 
example,  falls  into  the  above  class  if  orthogonality  is  not  enforced. 

A less  general  algebraic  transformation  (referred  to  as  an  independent 
algebraic  transformation)  composed  of  independent  transformations  in  each  direc- 
tion can  be  written  as 


0 g 


(2) 
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where  h and  g are  two  independent  parameters.  For  this  transformation,  dx  and 
dy  are  orthogonal  but  their  magnitudes  are  not  equal  when  the  magnitudes  of  d£ 
and  dn  are  equal . 

An  orthogonal  differential  coordinate  transformation  can  be  written  in  a 
similar  manner  as 


dx' 

h O' 

rr 

• 

dv 

.0  g. 

co  80 
sin0 


-sin0 

COS0 


d£ 

dn. 


(3) 


where  h,  g,  and  0 are  three  independent  parameters,  0 is  the  angle  between  the 
dx  and  d£  directions,  and  dx  is  orthogonal  to  dy. 

A conformal  differential  coordinate  transformation  can  be  written  as 


dx' 

dy 

where  h and  0 are  two  independent  parameters.  For  this  transformation,  dx  and 
dy  are  orthogonal  and  their  magnitudes  are  equal  when  the  magnitudes  of 
d£  and  dn  are  equal.  Thus  the  magnification  h (called  the  metric  or  mapping 
modulus)  between  the  (x,y)  and  (£,n)  planes  is  not  a function  of  direction;  it 
is  Just  a function  of  location.  A conformal  transformation  is  also  angle  pre- 
serving in  a local  sense,  so  that  if  two  differential  vectors  meet  at  the  angle 
a in  the  (£,n)  plane,  they  also  meet  at  the  angle  a in  the  (x,y)  plane. 


h 0 
0 h 


CO  80 

-sin© 

dS 

sinO 

COS0 

dn 

(4) 


FLOW  SIMULATION  IMPLICATIONS 

A grid  can  affect  the  storage  and  speed  of  a flow  simulation.  This  section 
touches  on  some  of  the  aspects  of  this  dependence  of  a flow  simulation  on  the 
grid. 

The  computer  storage  required  for  a grid  was  once  an  important  consideration 
for  flow  simulation.  For  a finite  difference  solution,  the  independent  para- 
meters in  the  differential  coordinate  transforms  above  were  usually  stored  (for 
an  orthogonal  grid)  or  regenerated  as  required  (for  a grid  constructed  using 
only  two  independent  algebraic  stretchings) . For  a limited  set  of  flow 
problems,  which  fortunately  includes  two-dimensional  transonic  potential  flow, 
it  turns  out  that  a conformal  mapping  requires  the  storage  of  only  the  metric, 
h,at  each  grid  point  in  addition  to  the  potential,  (f>.  The  x and  y physical 
plane  coordinates  and  the  parameter  0 in  Eq.  (4)  are  not  needed  throughout  the 
grid  in  this  case,  due  to  the  way  the  gradient  operator  transforms  under  a 
conformal  mapping,  so  only  two  quantities  need  be  stored  for  each  grid  point. 
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This  is  one  reason  why  conformal  grids  were  popular  in  the  early  days  of  com- 
puters for  finite  difference  solutions  using  a velocity  potential.  Purely 
algebraic  grids  were  also  popular  as  they  could  be  regenerated  as  required,  so 
that  only  the  potential  need  be  stored  at  each  grid  point.  The  rapidly  declin- 
ing cost  of  computer  memory  has  made  storage  requirements  a minor  issue*  at 
many  installations  for  two-dimensional  and  axisymmetric  inviscid  flow  problems; 
in  fact,  a modem  finite  volume  technique  may  routinely  store  twenty  to  thirty 
quantities  at  each  grid  point.  The  storage  requirements  of  a grid  system 
therefore  need  not  be  a deciding  factor  in  the  choice  of  a grid  system. 

The  computer  time  required  by  a flow  algorithm  depends  on  the  nature  of  the 
grid  system.  For  example,  when  a conformal  grid  (or  an  orthogonal  grid  using 
a conformal  step  followed  by  independent  algebraic  stretchings)  is  used,  a 
byproduct  is  the  solution  for  the  two-dimensional  incompressible  flow  as  dis- 
cussed later.  This  incompressible  solution  may  often  be  used  with  a compressi- 
bility transformation  as  an  accurate  first  guess  for  an  iterative  solution  of 
the  compressible  problem,  reducing  the  computer  time  needed  to  converge  the 
iterative  process.  In  addition,  a suitable  conformal  mapping  often  automati- 
cally provides  a two-dimensional  concentration  of  the  mesh  in  regions  of  high 
gradient  (because  the  mapping  modulus  is  independent  of  direction  for  a con- 
formal mapping),  without  simultaneously  producing  a one-dimensional  concentra- 
tion where  not  needed.  This  allows  the  use  of  fewer  grid  points,  which  results 
in  faster  calculations. 

DESIGN  IMPLICATIONS 

A number  of  modern  aerodynamic  design  methods.  Refs.  [2,3],  wherein  the 
pressure  is  specified  and  the  body  shape  giving  this  pressure  is  calculated, 
use  conformal  transformation  techniques.  The  computer  programs  required  for 
such  a conformal  mapping  based  aerodynamic  design  process  are  minor  variations 
of  the  conformal  mapping  programs  required  to  generate  an  orthogonal  grid. 

For  self-consistency,  one  might  as  well  use  the  same  routines  for  both.  This 
may  prove  to  be  a significant  justification  for  conformal  grid  generation. 

COMPLEX  VARIABLE  NOTATION 

Complex  variable  notation  allows  us  to  write  the  differential  conformal 
transformation,  Eq.  (4),  in  a compact  form  as 


* In  some  computing  installations,  the  charging  algorithm  recognizes  this  fact 
so  that  the  charge  for  a calculation  only  depends  on  the  CPU  time  and  not  on 
the  product  of  CPU  time  and  storage. 


Ill 


or 
or 
or 

where 
and 
and 

Thus  given  the  relation  Z = F(c)  we  can  easily  calculate  the  mapping  modulus  as 
dZ  f 

h = -j — = | F ' ( C ) | . Complex  notation  greatly  simplifies  the  use  of  conformal 
d C i 

transformation  techniques,  and  is  well  supported*  on  most  modem  computing 
systems. 

CONFORMAL  MAPPING 

Conformal  mapping  provides  a surface  correspondence  that  is  not  limited  to 
planar  surfaces.  Conformal  mappings  from  non-planar  surfaces,  such  as  fron 
spheres  or  general  cylindrical  surfaces,  to  a plane  have  been  employed  by 
cartographers  under  the  title  of  projections.  We  now  will  discuss  some  charac- 
teristics of  these  different  types  of  mappings,  starting  with  some  simple  planar 
mappings. 

Planar  mappings,  relating  points  in  two  different  planes,  are  the  most 
familiar.  An  important  planar  mapping  is  the  bilinear  transform  given  by 

Lz.Zg.)  (c_a)  = (c-A)  (5) 

(Z-b)  U a)  (C-B)  ^ 

where  the  points  Z * a,  b,  and  c correspond  to  the  points  £ « A,  B,  and  C 
respectively.  It  can  be  shown  that  Eq.  (5)  transforms  circles  or  straight  lines 
in  the  Z plane  to  circles  or  straight  lines  in  the  C plane.  A region  at 
infinity  can  be  treated  as  a point  in  this  equation.  This  mapping  illustrates 
that  one  may  specify  the  correspondence  of  only  three  specific  points  between 
planes  (i.e.,  one  may  assign  a,  b,  c and  A,  B,  C in  an  arbitrary  manner).  The 
bilinear  transformation  is  often  used  after  another  conformal  mapping  step  to 
produce  a canonical  (or  standardized)  contour  for  the  next  step. 


dx  + idy  = he^®(dC+idri) 
dZ  = H d? 

Z =/ H dC 
Z = F(C) 

Z = x + iy 
H = hei0 
C = 5 + in 


* A breakup  into  real  and  imaginary  parts,  and  then  using  real  arithmetic  only, 
can  produce  a two  to  one  speed  improvement  on  IBM  systems.  This  means  that 
inner  loops  for  heavily  used  routines  should  be  considered  for  such  a breakup. 
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Often  a contour  is  conformally  mapped  onto  a unit  circle  and  then  the  con- 
formal transformation 

u = £n(c) 


= Jln(re^) 


inr  + i0 


(6) 


is  applied, followed  by  the  independent  algebraic  transformation 

_ inr 
R - e = r 

0 = 6 


(7) 


to  map  the  interior  of  the  unit  circle  to  the  interior  of  the  rectangle  given 
by 

OtRsl 

0 £ 0 £ 2ir 


A uniform  grid  in  this  rectangle  plane  may  then  be  constructed  and  mapped  back 
to  the  physical  plane  where  the  resulting  grid  will  be  orthogonal. 

A simple  conformal  mapping  to  map  a general  axisymnetric  cylindrical  surface 
onto  a plane  has  been  variously  reported  in  the  literature  and  is  well  recog- 
nized overseas,  Refs.  [4,  5,  6,  7,  8,  9],  but  evidently  not  at  present  in  the 
United  States.  This  mapping  is  given  by 


dZ  = r d; 


(8) 


where 

dZ  = ds  + irdO 


dC  = d£  + idn 


and 


s = arc  length  along  cylinder  in  axial  direction 


i 

2 2 
] dx  . 


The  modulus  of  this  transformation  is  simply  r,  the  radius  in  the  cylindrical 
coordinate  system  (x,r,9)  describing  the  physical  surface.  For  a constant 
radius  axisymnetric  cylinder,  the  radius  r drops  out  of  the  flow  equations  and 
one  can  think  of  simply  unwrapping  the  cylinder,  but  in  fact  one  is  implicitly 
performing  a conformal  mapping.  For  an  axisymmetric  conical  surface,  the  above 
relation  becomes 

Z = eC8ln*  , (9) 
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where  $ is  the  angle  between  the  axis  and  the  cone  surface.  One  obvious  use 
of  cylindrical  mappings  is  in  the  generation  of  grids  for  quasi  three-dimensional 
blade-to-blade  flow  simulations  for  axial,  mixed  (meaning  both  significant 
radial  and  axial  flow  components),  and  centrifugal  flow  turbomachines.  For 
= it/2,  Eq.  (9)  reduces  to  Eq.  (6)  which  has  been  used  in  Ref.  [10]  to  map  a 
centrifugal  impeller  with  near  log-spiral  blades  onto  a two-dimensional  cascade 
of  blades. 

The  requirements  for  constructing  a useful  geographical  or  astronomical  map 
(i.e:  local  angles  are  preserved  and  the  local  scale  factor  does  not  depend 

on  orientation)  are  identical  to  the  definition  of  a conformal  mapping.  Thus, 
the  mappings  developed  by  cartographers  to  represent  a spherical  surface  by  a 
planar  map  apply  directly.  Among  such  maps  are  stereographic  and  Mercator 
projections,  and  a more  general  projection  used  for  star  maps,  Ref.  [9]. 
Stereographic  projections,  Invented  by  the  Greek  astronomer  Ptolemy  as  cited 
in  Ref.  [11],  were  used  to  develop  grid  systems  for  the  computation  of  super- 
sonic flow  over  conical  bodies.  Ref.  [12],  and  more  general  bodies.  Ref.  [13], 

A spiral  groove  spherical  bearing  lubrication  study.  Ref.  [14],  also  used 
stereographic  projection.  Stereographic  projections  have  been  known  for  13 
centuries;  the  recent  innovation  was  in  how  to  use  them  for  flow  calculations. 

RIEMANN  SHEET  DETERMINATION 

The  bane  of  conformal  mapping  is  the  possibility  of  getting  on  the  wrong 
Riemann  sheet,  or  equivalently,  choosing  an  inappropriate  root.  Development 
of  conformal  mapping  computer  programs  would  be  substantially  easier  were  it 
not  for  the  multivalued  nature  of  many  mappings.  This  section  discusses 
techniques  that  have  proven  useful  for  determining  the  appropriate  root. 

To  illustrate  the  problem,  consider  the  mapping 


where  k is  a real  number.  We  can  express  £ in  polar  coordinates  as 

r 1(0+2™) 


where  n is  any  integer,  since 


Then 


_ k ikQ  ik2wn 
Z « r e e 
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For  a square  root  transformation  k = 1/2,  and  the  e*-k2irn  term  is  1 if  n is  even 
and  -1  if  n is  odd  so  that  there  are  two  different  solutions  for  Z for  the  same 
C (while  there  is  only  one  value  of  C for  each  value  of  Z) . This  multivalued- 
ness can  be  visualized  in  terms  of  Riemann  sheets,  with  each  different  value  of 
C lying  on  a different  sheet  as  illustrated  in  Fig.  1.  For  k « 1/3,  there  are 
three  values  of  Z for  each  value  of  C.  For  irrational  values  of  k,  there  are 
an  infinite  number  of  values  of  Z for  each  value  of  c. 


Z = 


Z plane 


f plane 


FIGURE  1.  RIEMANN  SHEETS 


A computer  implementation  of  Z = q'  using  the  CSQRT  complex  square  root 
FORTRAN  function  will  return  only  the  root  having  a positive  real  part;  the 
other  root  will  be  the  negative  of  this  root.  It  is  up  to  the  investigator  to 
provide  program  logic  to  choose  which  of  these  two  roots  is  appropriate. 

For  k *=  1/2,  which  is  rather  common,  physical  reasoning  can  often  be  used 
to  select  the  appropriate  value  of  Z.  In  particular,  mappings  using  the 
principle  of  reflection  (covered  later)  produce  quadratic  equations  with  two 
roots  which  are  image  points  with  respect  to  a circle.  In  this  case  one  may 
choose  the  root  either  inside  or  outside  the  circle  based  on  physical  reasoning. 
For  the  square  root  transformation,  as  often  used  in  airfoil  or  wing  grid 
generation,  or  the  transformation, 


Z = £ + e 


(11) 
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which  also  has  two  roots  and  is  often  used  for  inlet  grid  generation,  the  root 
selection  is  rather  simple.  Referring  to  Figs.  2 and  3,  the  initial  value  at 
"A"  is  the  root  with  a negative  real  value,  and  as  one  traverses  the  airfoil  in 
a clockwise  direction  and  the  inlet  in  a counterclockwise  direction,  the  root 
to  the  right  of,  and  (of  those  remaining,  if  any)  nearest  to,  the  previous  root 
is  chosen.  This  simple  algorithm  has  worked  well  for  a large  variety  of  airfoil 
and  inlet  grid  calculations. 


Z plane  Z plane 


t plane 


B 


_rr^p_ 

FIGURE  2.  AIRFOIL  MAPPING 


f plane 


FIGURE  3.  NACELLE  MAPPING 


For  a transformation  involving  an  irrational  power  of  k near  (but  not 
exactly)  1/2,  there  are  many  values  of  Z that  lie  near  each  other  and  the 
selection  is  less  obvious.  The  von  Kartnan-Treff tz  transformation,  Ref.  [15], 
given  by 


which  transforms  an  airfoil  in  the  Z plane  into  a near  circle  in  the  C plane, 
is  one  such  case.  When  k *■  ^ , the  trailing  edge  is  cusped  and  Eq.  (12) 
becomes  the  Joukowski  transformation. 

One  technique  to  deal  with  the  root  problem  is  that  of  "tracking",  where  one 
chooses  the  value  of  n such  that  (0+2Trn)  in  Eq.  (10b)  is  closest  to  the  value 
of  (0+27rn)  for  a nearby  point.  The  value  of  0 is  obtained  from  the  real  and 
imaginary  parts  of  C using  the  FORTRAN  ATAN2  function,  which  returns  a value 
of  0 between  -it  and  tt.  In  this  manner  one  "tracks"  the  argument  (9+2un)  con- 
tinuously along  a curve  and  chooses  n to  make  this  argument  continuous.  There 
remains  the  problem  of  how  to  choose  the  root  for  the  initial  point.  This  may 
often  be  accomplished  from  knowledge  of  the  behavior  of  the  transformation  at 
the  point  in  the  mapped  plane  which  corresponds  to  the  far  field  in  the  physical 
plane. 
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Another  method  of  root  selection  is  to  introduce  a "branch  cut"  based  on 
physical  reasoning,  and  then  preassign  n on  each  side  of  the  cut.  The  choice 
of  root  for  a reflection  mapping  as  described  earlier  is  equivalent  to  using  a 
circle  for  the  cut.  For  many  problems  the  cut  may  be  taken  as  a straight  line 
and  is  chosen  so  as  to  not  cross  the  boundary  being  mapped.  Some  mappings  may 
require  more  than  one  cut,  as  illustrated  in  Refs.  [16]  and  [17].  Using  a 
branch  cut  can  be  considered  as  a static  technique  to  preassign  n on  a global 
basis,  while  tracking  is  a dynamic  technique  to  determine  n on  a local  basis. 

It  is  not  unusual,  for  complicated  mappings,  to  occasionally  encounter  root 
selection  errors  in  a conformal  mapping  program  which  had  been  thought  to  be 
resistant  to  such  effects.  This  usually  happens  when  attempting  to  map  a 
geometry  significantly  different  from,  but  still  in  the  same  class  as,  the 
geometry  used  to  develop  the  program.  For  simple  mappings  (i.e:  those  quoted 

so  far  in  this  discussion  with  the  possible  exception  of  Eq.  (12))  root  selec- 
tJrn  problems  ere  not  expected  if  the  above  techniques  have  been  carefully 
implemented.  In  any  case,  an  error  is  easily  detected  by  defining  a uniform 
grid  in  the  rectangular  computational  plane  and  mapping  this  grid  back  to  the 
physical  plane  where  it  is  plotted.  If  the  physical  plane  grid  is  continuous, 
does  not  overlap,  and  does  not  exhibit  open  sections  or  jumps,  the  roots  have 
probably  been  selected  correctly.  It  is  important  to  visually  inspect  the  grid 
created  for  each  new  geometry  before  it  is  used  for  a flow  calculation  if  there 
is  any  doubt  on  this  point. 

CLASSICAL  MAPPING 

The  classical  technique  of  mapping  consists  of  mapping  a contour  (an  airfoil, 
a cascade,  an  inlet,  etc.)  to  a near-circle  by  a single  transformation  or  a 
sequence  of  simple  transformations.  The  near-circle  is  then  mapped  to  a circle 
by  a near-circle  to  circle  transformation,  such  as  proposed  by  Theodorsen  and 
Garrick,  Ref.  [18]. 

The  curve  to  near-circle  transformation(s)  often  require  a substantial  level 
of  ingenuity  and  luck,  and  even  then  can  fall  prey  to  root  selection  errors. 

In  simple  cases,  such  as  for  an  isolated  airfoil  or  an  inlet,  the  near-circle 
is  actually  near  to  being  a circle.  For  other  cases  (inlet  with  centerbody  far 
inside  or  far  outside,  turbine  cascade,  small  gap/chord  compressor  cascade),  the 
near-circle  may  not  even  resemble  a circle.  When  this  happens,  the  investigator 
must  either  find  a better  mapping  to  a near-circle  (which  may  not  be  easy  in 
practice,  even  after  many  years  of  experience),  or  change  to  a nonclassical 
technique  such  as  discussed  later  in  the  section  on  one-step  mappings. 
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SEQUENTIAL  MAPPINGS 

For  a complicated  mapping  of  the  form  Z - f(c).  roots  are  often  more  easily 
selected  and  the  mapping  Is  more  easily  created.  If  the  process  Is  broken  up 
Into  an  equivalent  sequence  of  simpler  mappings.  For  Instance,  the  von  Karman- 
Trefftr  transformation  of  Eq.  (12)  can  be  restated  as  the  sequence 


u 

n 


Z-A 

Z-B  ’ 


C 


(13a) 

(13b) 

(13c) 


The  bilinear  transformations,  Eqs.  (13a)  and  (13c)  require  no  root  selection, 
while  Eq.  (13b)  requires  keeping  track  of  only  a single  angle.  Other  schemes 
may  require  the  separate  tracking  of  the  angles  related  to  both  the  (Z-A)  and 

k 

the  (Z-B)  term. 

The  logarithmic  transformation,  Eq.  (6),  and  Its  Inverse 


C ■ eZ  (14) 

are  the  key  elements  of  conformal  mappings  for  a cascade  of  airfoils,  such  as 
occur  In  turboanchlnery.  For  this  case,  there  are  an  Infinite  number  of 
Identical  airfoils  in  the  Z plane,  all  of  which  map  onto  the  same  contour  In 
the  t plane.  When  mapping  from  the  c plane  to  the  Z plane,  the  correct  root 
Is  the  one  that  "tracks"  a nearby  point;  the  other  roots  will  be  different  by 
an  Integral  multiple  of  2*  In  the  vertical  (or  Imaginary  axis)  direction.  Most 
of  the  cascade  mappings  which  the  author  has  seen.  Refs.  [19] . [20],  [21]  and 
[22],  can  be  restated  as  a simple  sequence*  with  Eq.  (14)  (possibly  combined 
with  a bilinear  transformation)  as  the  first  element. 

As  a general  principle,  It  Is  highly  recomsended  for  mappings  of  the  form 
Z - f(c)  that  a complicated  mapping  be  restated  as  a sequence  of  simpler  map- 
pings for  actual  implementation  in  a computer  program. 


NEAR-CIRCLE  TO  CIRCLE  MAPPING 

If  an  orthogonal  grid  Is  desired,  the  mapping  of  a near-circle  to  a circle 
Is  often  chosen  as  one-step  of  the  mapping  process.  Although  there  are  now  a 
plethora  of  ways  to  accomplish  this  mapping**,  the  classical  technique  is  that 
of  using  the  Theodorsen-Garrlck  transformation.  Ref.  [18],  given  by 

* The  author  has  not  yet  been  successful  In  breaking  up  the  Garrick  transform. 
Ref.  [23],  Into  a simple  form. 

**  Including  panel  techniques,  Schwartz-Chrlstoffel  techniques,  and  elliptic 
techniques  with  orthogonal  boundary  control. 
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j=N 

£(a  +ib  Kj 
Z 1=0  3 : 

c *e 

or,  in  its  derivative  form  by 

j=N 

^(a,+ib 
dZ  j =0  J J 
d?  = 6 


(15a) 


(15b) 


What  is  frequently  not  appreciated  is  that  when  Eqs.  (15a,  15b)  are  implemented 
in  the  required  iterative  manner  (see  Ref.  [18  or  24]  for  details),  the  itera- 
tion may  not  necessarily  converge  unless  the  near-circle  is  sufficiently  "near" 
to  being  a circle.  In  particular,  there  is  a requirement  on  the  maximum  value* 


allowed  for 


d(inr) 


d0 


for  the  mapping  in  Eq.  (15a),  where  r and  0 are  the  polar 
coordinates  describing  the  near-circle.  For  the  vast  majority  of  airfoils,  the 
classical  von  Karman-Tref ftz  transformation  produces  a near-circle  that  meets 
the  above  requirement.  For  turbine  cascades,  on  the  other  hand,  the  above 
requirement  is  usually  not  met.  It  is  possible  to  underrelax  the  iteration, 
but  the  required  underrelaxation  parameter  may  need  to  be  quite  small  for 
Eq.  (15a)  to  ensure  convergence,  resulting  in  unrealistic  computation  time 
requirements.  The  form  given  by  Eq.  (15b)  is  much  less  sensitive,  and  with  an 
underrelaxation  factor  of  .5  has  converged  even  for  a case  where  r was  a 
multi  pie- valued  function  of  0,  so  that  was  infinite!  Obviously,  the 

convergence  criterion  for  Eq.  (15b)  is  substantially  different  than  that 
derived  for  Eq.  (15a)  in  Ref.  [25].  In  short.,  if  one  uses  the  derivative  form, 
the  chances  for  success  are  much  higher.  Of  course,  the  derivative  form 
requires  the  integration  of  c complex  function,  but  an  excellent  technique 
using  a spline  exists  as  discussed  later.  Both  of  the  above  forms  can  employ 
fast  Fourier  techniques,  as  in  Refs.  [24,  26],  to  efficiently  evaluate  the 
series  JZ(aj+i.bj ) at  evenly  spaced  increments  on  the  unit  circle  in  the  C 
plane  and  thus  keep  the  required  computing  time  modest. 


ONE-STEP  MAPPINGS 

A one-step  mapping  maps  a contour  onto  a canonical  shape,  such  as  a circle 

or  the  real  axis,  in  a single  step.  Such  transformations  are  usually  given  in 

dZ 

derivative  form,  namely  — * fU). 

d; 


i 


* For  a relaxation  parameter  of  unity,  this  maximum  value  is  .2955. 
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The  classic  Schwarz-Christof fel  transformation  for  a polygon  with  H straight 
sides  takes  the  form 


dZ 

dC 


j=N  k , 

T7  <?-b  ) 1 

j=l  J 


(1C) 


while  for  a polygon  with  curved  sides,  Davis,  Ref.  [27],  uses  a differential 
form  of  the  product  term  in  Eq.  (16)  so  that 


dZ  */^-b)dB 
dC  ' 6 


(17) 


where  B is  an  angle  related  variable.  Davis  employs  a composite  integration 
formula  to  resolve  the  curvature  effects.  Equation  (16),  as  used  by  Skulsky, 
Ref.  [28],  can  be  considered  a subset  of  the  Davis  technique.  The  Davis  tech- 
nique has  been  successfully  employed  in  Refs.  [12],  [29],  and  [30]  to  map  a 
wide  range  of  complicated  configurations. 

An  alternative  is  to  use  the  form 


j=N 

V 


dZ 

dC 


g(C)e 


j^<aj+lbi)s 


j 


(18) 


where  g(t)  is  chosen  to  resolve  angles  or  general  behavior*,  while  the  expo- 
nential term  accounts  for  the  curvature.  This  form  can  take  advantage  of  fast 
Fourier  techniques  to  evaluate  the  coefficients  a^  and  b^. 

One  example  of  such  a mapping  is  given  by 


dZ 

d£ 


(19) 


which  is  used  in  Ref.  [26]  to  map  an  airfoil  in  the  Z plane  to  a circle  in  the 
C plane.  Another  mapping  of  this  type  is  given  by 


* The  g(c)  function  is  easi’  v constructed  using  the  Schwarz-Christof  fel  tech- 
nique, as  illustrated  Ln  Kober,  Nehari,  or  Milne -Thomson,  Refs.  [31],  [32] 
or  [33]. 
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where 


and 


(;-iso)  Sw1 

2 eJ 

(I-IRqHC+I)  (c-l) 


Rq  = 1 + 2^te 

nose  radius 

0 X • - - ii... 

Inlet  interior  radius 


(20) 


This  mapping  was  developed  by  the  author  to  map  the  region  exterior  to  a semi- 
infinite  inlet  and  above  the  centerline  to  the  interior  of  a circle  as  illu- 
strated in  Fig.  4.  A similar  form,  suggested  by  Jameson,  Ref.  [34],  given  by 


dZ 

dC 


(4+1)  3(1-?) 


(21) 


was  used  to  map  an  inlet  and  centerline  to  a half  circle  as  Illustrated  in 
Fig.  5. 


tplWM 


FIGURE  4.  INLET  TO  CIRCLE  MAPPING  FIGURE  5. 

INLET  TO  HALF  CIRCLE  MAPPING 


The  latter  mapping  did  not  adequately  resolve  the  exterior  inlet  contour  far 
downstream  (but  was  good  otherwise),  while  a mapping  based  on  Eq.  (20)  ade- 
quately resolved  the  inlet  in  all  areas.  The  reason  is  that  the  evenly  spaced 
point  distribution  (as  required  by  fast  Fourier  techniques)  on  the  4 plane 
circle  is  too  sparse  when  transformed  back  to  the  Z plane  to  accurately  repre- 
sent the  nacelle  far  downstream  for  Eq.  (21).  This  example  illustrates  the 
need  to  choose  the  correct  canonical  domain.  Of  course,  the  two  canonical 
domains  are  related  in  a simple  quadratic  manner  through  the  transformation 


Z-l 

Z+l 


2 


(22) 
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which  transforms  a unit  circle  in  the  Z 
to  the  real  axis  from  -1.  to  1.)  in  the 

Z plan* 


plane  to  a unit  upper  half  circle  (and 
C plane  as  illustrated  in  Fig.  6. 

f plane 


£ 6. 

CIRCLE  MAPPING 


The  one-step  mapping  of  Eq.  (18)  is  far  simpler  to  program*  than  the  con- 
ventional classic  technique  described  earlier,  converges  stably  and  rapidly, 
and  is  easy  to  modify  later  for  a new  class  of  geometries.  If  one  wishes  to 
use  fast  Fourier  techniques  for  conformal  mapping,  Eq.  (18)  is  recommended  as 
the  method  of  choice  rather  than  using  a sequence  resulting  in  a near-circle 
followed  by  a derivative  transform  like  Eq.  (15b),  especially  since  Eqs.  (18) 
and  (15b)  differ  only  by  a very  few  lines  of  computer  code. 

It  is  reconmended  that  one  should  not  expect  the  Fourier  series  in  Eqs.  (15a), 
(15b),  or  (18)  to  resolve  a slope  discontinuity  by  the  sheer  brute  force  of  a 
large  number  of  Fourier  terms.  One  might  "get  away"  with  doing  so  for  small 
discontinuities  in  slope  (e.g.,  5 degrees)  where  a high  accuracy  near  the  dis- 
continuity is  not  required,  but  for  larger  slope  discontinuities  the  accuracy 
and  convergence  of  the  transformation  will  suffer.  It  is  easy  to  remove  such 
slope  discontinuities  using  the  hinge  point  transformations  covered  next  or  by 
incorporating  slope  discontinuities  into  the  g(t)  term  in  Eq.  (18)  so  that  the 
Fourier  series  only  has  to  resolve  a function  with  a continuous  derivative. 

HINGE  POINT  TRANSFORMATIONS 

At  the  other  extreme  from  one-step  mappings,  the  hinge  point  transformations 
due  to  Moretti  and  Hall,  Refs.  [17]  and  [35],  break  the  problem  up  into  a 
large  number  of  sequential  applications  of  a single  mapping.  For  problems  with 
a plane  of  symmetry,  Moretti  applied  a von  Karman-Tref ftz  transformation  to  each 
pair  of  slope  discontinuities  in  turn,  to  produce  a smooth  near-circle.  Hall 


* Note  that  the  computer  program  in  Ref.  [26]  can  be  easily  modified  to  imple- 
ment a mapping  using  Eq.  (18),  instead  of  using  Eq.  (19)  for  the  airfoil. 
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repeatedly  applied  the  transformation  Z = (C~Cq)  (Eq.  (10a)),  where  Cq  and  k 
have  been  chosen  to  remove  the  left-most  remaining  slope  discontinuity  at  each 
stage,  to  produce  a smooth  near  half  plane.  Once  a near-circle  or  near  half 
plane  is  obtained,  a non-conformal  shearing  or  a conformal  mapping  can  be  used 
to  produce  a canonical  domain. 


MULTIPLE  BODIES 

Techniques  to  map  multiple  bodies  can  be  classified  as  simultaneous,  sequen- 
tial, iterative,  and  periodic.  This  section  discusses  these  techniques. 

A simultaneous  multiple  body  mapping  maps  two  or  more  bodies  simultaneously 
to  near-canonical  or  canonical  domains.  The  single  mapping 


j=N 

n ^ 

J-l  1 Nj 


£=N 

TT 

£=1 


z~hi] 

Z"ZJ 


(23) 


taken  from  Ref.  [24]  simultaneously  maps  N airfoils  to  N near-circles.  The 
Garrick  transformation  in  Ref.  [36]  simultaneously  maps  two  concentric  near- 
circles to  two  concentric  circles. 

A sequential  multiple  body  mapping  maps  two  contours  to  canonical  contours 
by  first  mapping  one  contour  to  a canonical  contour,  then  mapping  the  second 
contour  to  a canonical  contour  while  preserving  the  nature  of  the  first 
canonical  contour.  Two  examples  of  this  are  given  in  Ref,  [24],  where  trans- 
formations are  given  to  map  an  airfoil  to  a near-circle  (or  to  map  a near-circle 
to  a circle)  while  keeping  a nearby  circle  a circle.  Another  example  is  illu- 
strated in  Fig.  7,  where  a centerbody  is  mapped  onto  the  real  axis  using  a 
technique  similar  to  that  in  Ref.  [37],  and  then  a nearby  inlet  is  mapped  to  a 
circle,  using  the  Eq.  (20),  while  keeping  the  real  axis  a straight  line. 


FIGURE  7.  INLET/ CENTERBODY  TO 
CIRCLE  MAPPING  SEQUENCE 
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After  doing  a number  of  problems  both  simultaneously  and  sequentially,  the 
author  has  found  that  a sequential  multiple  body  technique  usually  works  better 
in  practice  than  a simultaneous  technique  since  the  branch  cuts  required  to 
resolve  the  root  selection  problem  are  simpler.  This  more  than  balances  the 
increase  in  computer  time*  required  for  a sequential  approach. 

An  iterative  technique  to  map  multiple  contours  to  canonical  contours  has 
recently  been  developed  by  Halsey,  Ref.  [38],  In  this  technique  each  contour 
is  individually  mapped  to  a circle  in  sequence,  with  no  special  constraints 
on  holding  the  other  contour  shapes.  By  repeatedly  cycling  through  all  the 
contours  one  at  a time,  the  contours  are  all  mapped  to  circles  to  engineering 
accuracy.  That  this  process  converges  is  not  particularly  surprising;  what 
is  surprising  is  that  only  about  five  cycles  through  all  the  contours  are 
required  to  achieve  four  digit  accuracy  for  practical  cases. 

A periodic  mapping  maps  a finite  or  infinite  number  of  identical  (but  other- 
wise displaced  or  rotated)  contours  onto  overlayed  (but  otherwise  identical) 
contours  on  different  Riemann  sheets.  As  an  example,  the  logarithmic  trans- 
formation of  Eq.  (6),  when  appropriately  scaled,  maps  a cascade  of  airfoils  in 
the  u plane  onto  a single  (highly  distorted)  airfoil-like  contour  in  the  C 
plane.  The  transformation  Z « (which  is  a subset  of  Eq.  (10a)),  where  N is 

an  integer,  maps  a region  with  N angular  periodic  boundaries  in  the  Z plane  onto 
a region  with  one  periodic  boundary  in  the  C plane.  By  these  means,  a periodic 
configuration  can  be  reduced  to  a form  wherein  the  periodicity  condition  reduces 
essentially  to  a continuity  condition  (except  for  circulation-type  terms 
associated  with  branch  cuts). 

CREATING  NEW  MAPPINGS 

This  section  includes  a discussion  of  a few  general  techniques  that  can  be 
used  to  create  new  conformal  mappings,  and  some  restrictions  to  keep  in  mind 
while  doing  so. 

The  first  step  when  considering  the  development  of  a new  mapping  is  to 
verify  that  the  desired  mapping  has  not  already  been  created.  The  book  by 
Kober,  Ref.  [31],  includes  an  extensive  survey  of  mappings  developed  up  to  the 
mid  1940's,  and  should  be  reviewed  first.  The  reference  material  in  Table  1, 
included  later  in  this  paper,  can  be  consulted  by  topic.  If  none  of  the  above 
include  the  desired  mapping,  a careful  literature  search  may  be  warranted.  Only 
when  the  investigator  is  reasonably  sure  of  not  "reinventing  the  wheel"  should 
a new  mapping  be  attempted. 

* A 2-body  sequential  approach  may  typically  require  twice  the  computer  time  of 
a simultaneous  approach. 
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If  the  complex  potential  u for  the  two-dimensional  incompressible  flow  over 
(or  through)  a body  is  known,  i.e.,  if 

u - 4>  + i<|i  - f(Z) 

is  known,  then  the  mapping  from  the  Z plane  to  the  m plane  (with  coordinates 
(♦,'('))  is  known.  Thus, knowing  the  incompressible  flow  is  equivalent  to  knowing 
the  mapping  from  the  Z plane  to  a rectangle  in  the  u plane,  and  vice  versa. 

This  principle  can  be  used  to  construct  conformal  mappings.  In  particular,  a 
"panel  method"  potential  flow  program  can  be  used  to  construct  a conformal 
mapping,  Ref.  [39]. 

Another  technique  for  the  creation  of  new  mappings  can  be  summarized  by  the 
description  "guess  and  plot."  Using  this  technique  a function  having  appro- 
priate zeroes,  poles,  and  singularities  is  guessed*,  and  its  effect  is  deter- 
mined by  plotting  a transformed  contour  or  a transformed  grid.  About  one  third 
of  the  functions  this  author  has  guessed  over  the  years  have  had  the  appropriate 
action,  indicating  that  this  technique  is  more  viable  than  might  be  thought  at 
first  glance. 

Since  a conformal  mapping  is  simply  a functional  relationship,  if  one  knows 

Z - f(w)  , 

and 

C - g(“) 

then  one  has 

z - f(g-1(0)  . 

In  short,  if  one  knows  how  to  map  two  different  contours  to  the  same  contour, 
then  one  can  map  one  of  these  contours  onto  the  other,  Refs.  [3i]  and  [40], 

Put  another  way,  a conformal  mapping  of  a conformal  mappirp  is  a conformal 
mapping.  One  then  begins  to  think  in  terms  of  building  up  a sequence  of  map- 
pings, with  each  step  bringing  a contour  closer  to  a canonical  contour  such  as 
a unit  circle.  Thus  by  combining  known  mappings  in  a new  sequence,  one  can 
construct  new  mappings. 

In  generating  an  orthogonal  grid  with  conformal  mapping,  there  are  two 
processes  involved.  The  first  process  is  to  conformally  map  the  desired  con- 
tour (i.e:  an  airfoil,  or  inlet,  etc.  ...)  to  a rectangle  (or  equivalently 

to  a circle  or  half  plane).  The  computing  time  required  for  this  step  is  often 

* These  singularities,  poles,  and  zeroes  are  placed  at  points  on  the  contour 
where  the  slope  is  discontinuous,  and  when  necessary  (e.g.,  near  regions  of 
high  convex  curvature  of  the  contour),  inside  the  contour  along  a line  join- 
ing the  local  contour  point  of  maximum  curvature  and  its  center  of  curvature. 

How  far  along  this  line  depends  on  the  total  slope  change  nearby,  and  the 
singularity  powers  also  depend  on  these  slope  changes. 
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nearly  linearly*  proportional  to  L,  the  number  of  input  points.  The  second 
process  is  classically  to  map  an  orthogonal  rectangular  grid  in  the  mapped 
plane  back  to  the  physical  plane.  This  involves  a computing  time  proportional 
to  MxN,  the  total  number  of  grid  points.  It  is  common  for  MxN  to  be  substan- 
tially larger  than  L.  This  implies  that  the  forward  transformations  can  be 
somewhat  inefficient  (l.e.,  implicit)  but  the  inverse  transformations  should 
be  efficient  (i.e.,  explicit),  if  possible.  This  should  be  kept  in  mind  when 
creating  mappings  that  are  to  be  used  to  transform  a grid  back  to  the  physical 
plane. 

When  creating  new  mappings,  one  principle  to  be  followed  is  "do  not  create 
singularities**  within  the  flow  field";  however,  singularities  are  often  neces- 
sary and  completely  acceptable  within  a body,  on  a bounding  surface,  or  at 
images  of  "infinity."  An  example  of  creating  a singularity  within  the  flow 
field  is  shown  in  Fig.  8,  taken  from  Ref.  [19].  Figure  9 from  this  same 
reference  shows  a grid  system  with  the  same  number  of  grid  points  which  is 
better  suited  for  flow  computations.  The  above  principle  was  used  In  Ref.  [24] 
to  determine  some  mapping  parameters  that  were  otherwise  undetermined.  With 
another  choice  of  parameters,  the  mappings  therein  looked  something  like  those 
in  Fig.  10  rather  than  like  those  in  Fig.  11  where  the  parameters  were  chosen 
to  obey  this  principle. 


FIGURE  8.  GRID  SYSTEM  WITH  FIGURE  9.  GRID  SYSTEM  WITHOUT 

SINGULARITY  IN  FLOW  FIELD  SINGULARITY  IN  FLOW  FIELD 


* Actually  the  operation  count  goes  as  Lxl^L,  as  dictated  by  using  fast  Fourier 
techniques  for  Eqs.  (15a),  (15b),  or  (18). 

**  A singularity  occurs  when  kk  = 0 or  infinity. 
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Z plane 


f plane 


Z plane  t plan* 


FIGURE  10.  INCORRECT  CHOICE  OF 
MAPPING  CONSTANTS 


FIGURE  11.  CORRECT  CHOICE  OF 
MAPPING  CONSTANTS 


REFLECTION  PRINCIPLE 

It  is  possible  to  map  a contour  to  a desired  shape  while  simultaneously 
keeping  a nearby  straight  line  straight  by  using  a reflection  principle.  If 
the  relation 


f(Z)  = g(c) 

conformally  maps  a contour  from  the  Z plane  to  the  C plane  in  a desirable 
manner,  then  the  relations 


f(Z)-f(Z)  - g(C)-g(c)  , 

g(?) 


(24a) 

(24b) 


where  the  superscript  denotes  the  complex  conjugate  operator,  will  both 
have  a "similar"*  effect  but  will  preserve  the  shape  of  the  real  axis.  A 
proof  of  the  real  axis  shape  preservation  is  outlined  in  Appendix  A.  Examples 
involving  mappings  symmetric  with  respect  to  a line  appear  in  Refs.  [16]  and 
[17].  An  extension  of  the  above  concept  allows  the  construction  of  mappings 
which  preserve  a circle  by  utilizing  operators  involving  image  points  with 
respect  to  circles  rather  than  the  complex  conjugate  operator.  Examples  of 
circle  preserving  mappings  are  contained  in  Refs.  [24]  and  [36].  The  question 
of  whether  to  use  the  product  mapping,  Eq.  (24a),  or  the  ratio  mapping, 

Eq.  (24b),  can  be  resolved  by  simply  programming  both  and  then  choosing  the 
one  which  produces  the  most  desirable  effect. 


* "Similar"  is  taken  here  to  mean  not  identical,  but  generally  of  the  same 
nature. 
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CALCULATING  THE  GRID 

Recently,  Sockol  and  Adamczyk  in  Refs.  [39]  and  [41]  have  observed  that  the 
generation  of  an  orthogonal  grid  can  be  broken  down  into  two  independent  steps, 
namely; 

1)  calculation  of  the  boundary  point  correspondence  between  the 
physical  and  computational  planes  by  means  of  conformal  mapping 
or  some  other  orthogonality  preserving  technique,  and 

2)  the  generation  of  a grid  given  this  boundary  point  correspondence. 

This  observation  is  to  conformal  grid  generation  as  the  finite  volume  (or 

finite  element)  technique  is  to  flow  calculations;  it  breaks  the  grid  genera- 
tion problem  up  into  two  nearby  independent  steps  just  as  the  finite  volume 
technique  breaks  a flow  simulation  up  into  two  nearby  independent  steps 
(namely  the  generation  of  a grid  and  the  solution  of  the  flow  equations  on  the 
grid).  This  allows  a great  deal  of  flexibility  to  be  attained.  In  short,  what 
is  the  most  efficient,  simple,  or  general  technique  required  to  complete 
step  #1  (which  is  basically  a boundary  operation)  is  not  necessarily  the  same 
as  the  most  efficient,  simple,  or  general  technique  required  to  complete 
step  If 2 (which  is  basically  an  area  operation).  This  represents  a major  step 
forward  in  orthogonal  grid  generation.  This  section  reviews  some  methods  used 
to  complete,  step  It2. 

One  obvious  way  to  generate  an  orthogonal  grid  is  to  invert  the  conformal 
mapping  used  in  step  If  1,  and  then  map  an  orthogonal  grid  constructed  in  the 
computational  plane  back  to  the  physical  plane,  as  mentioned  earlier.  If  one 
has  used  an  algebraic  stretching  in  the  direction  tangential  to  the  body  to 
create  the  computational  grid,  fast  Fourier  techniques  cannot  be  easily  employed 
in  the  mapping  inversion.  Depending  on  the  mapping  complexity,  the  number  of 
Fourier  terms,  and  the  presence  of  implicit  expressions  in  the  inverse  mapping, 
the  creation  of  a grid  at  NxM  points  can  require  the  order  of  CxLxNxM  operations, 
where  L is  the  number  of  Fourier  terms  and  C is  a relatively  large  constant. 

In  practical  terms,  this  can  require  about  six  seconds  of  time  on  an  IBM  3081 
computer  for  a typical*  case.  For  a non-stretched  isolated  airfoil  mapping, 
where  fast  Fourier  techniques  can  be  used  in  the  inverse  mapping,  a typical 
grid  can  be  generated  in  about  two  seconds  of  IBM  3081  computer  time,  once  the 
conformal  mapping  is  known. 


* For  the  purposes  of  this  paper  "typical"  means  a 128x32  grid  using  256  Fourier 
terms  for  an  inlet/centerbody  configuration,  unless  otherwise  specified.  The 
IBM  3081  computing  time  quoted  can  be  taken  as  roughly  equivalent  to  an  equal 
amount  of  CDC  7600  computer  time. 
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Another  way  to  generate  an  orthogonal  grid  is  to  solve  a Laplace  problem 
(with  Dirichlet  boundary  conditions)  on  the  computational  grid  for  both  x and  y 
(the  physical  plane  coordinates),  which  are  known  on  the  boundary  as  described 
in  Refs.  [39]  and  [41].  Since  the  Laplacian  operator  will  retain  orthogonality, 
the  (x,y)  values  thus  calculated  will  describe  an  orthogonal  grid  in  the 
physical  plane.  A fast  Poisson  solver.  Ref.  [42],  with  an  operation  count  of 
CxNxMxS^N,  where  C is  a relatively  small  constant,  may  be  used  to  solve  the 
Laplace  problem.  For  a typical  grid,  this  will  require  about  two  seconds  on  an 
IBM  3081  computer.  This  solver  is  sufficiently  general  even  for  use  when 
independent  stretchings  are  used  following  a conformal  mapping  to  generate  the 
rectangular  grid  in  the  computational  plane. 

THREE  DIMENSIONS 

Conformal  mapping  is  basically  a surface  technique,  but  it  can  be  used  as 
one  component  of  a three-dimensional  grid  generation  system.  One  example  is 
shown  in  Fig.  12,  taken  from  Ref.  [16],  where  a three-dimensional  inlet/ 
centerbody  grid  was  constructed  using  conformal  mapping  in  each  circumferential 
plane,  or  slice,  followed  by  independent  algebraic  stretchings  to  construct  an 
orthogonal  grid  in  the  sliced  plane,  one  of  which  is  shown.  The  resultant 
three-dimensional  grid  is  orthogonal  in  two  directions,  but  not  in  the  third. 
Similar  techniques  are  used  to  construct  near-orthogonal  grids  in  Refs.  [13], 
[43],  and  [44]. 

REFERENCE  MATERIAL 

The  following  table  summarizes  a limited  amount  of  reference  material  for 
conformal  mapping  of  aerodynamic  configurations.  The  book  by  Kober,  Ref.  [31], 
is  particularly  useful  when  attempting  to  map  a new  configuration.  The  paper 
by  Moretti,  Ref.  [45],  is  highly  recommended  as  an  alternate  review  of  the 
field. 


FIGURE  12.  PERSPECTIVE  VIEW  OF 
THREE-DIMENSIONAL  COORDINATE  SYSTEM 


TABLE  1 - A LIMITED  SET  OF  CONFORMAL  MAPPING  REFERENCES 


Sublect 

References 

General  - Papers 

27,  28,  40,  45,  46,  47,  48 

General  - Books 

11,  31,  32,  33 

Airfoil 

18,  24,  26 

Airfoil  with  Spoiler 

49 

Two  Piece  Airfoil 

24,  36,  38 

Multiple  Bodies 

24,  38 

Inlet 

16,  50,  51 

Inlet  with  Centerbody 

16 

Axisymmetric  Body 

37 

Cascade 

19,  20,  21,  22,  23,  39,  41 

Non- Planar 

4,  5,  6,  7,  8,  9,  10,  12,  13,  14 

Symmetric 

16,  17 

Ducts 

27,  29,  30 
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PACKAGED  TOOLS 

To  accomplish  simple  conformal  transformations,  one  does  not  need  much  in 
the  way  of  supporting  tools.  For  more  complex  transformations,  a number  of 
packaged  tools  can  substantially  simplify  program  development.  Among  these 
tools  are  spline  fits,  fast  Fourier  transforms,  efficient  matrix  techniques, 
fast  Poisson  solvers,  near-circle  to  circle  transformation  packages,  and  good 
graphics.  The  need  for,  use,  and  availability  of  these  packages  is  covered  in 
this  section. 

Mapping  procedures  often  require  interpolation  between  points.  A spline  fit 
is  a good  way  to  accomplish  such  an  interpolation.  Ref.  [52].  In  addition,  it 
may  be  necessary  to  numerically  integrate  a function.  Such  an  integration  can 
be  accomplished  by  analytic  integration  of  a spline  passed  through  the  function 
values.  A spline  routine  accomplishing  both  of  these  objectives  is  listed  in 
a FORTRAN  program  form  on  pages  279-281  of  Ref.  [26].  By  simply  declaring  all 
floating  point  variables  in  this  routine  as  COMPLEX,  this  program  is  suitable 
for  integrating  a complex  function.  The  derivative  form  of  transformations  of 
Eqs.  (15b),  (16),  (17),  and  (18)  require  such  an  integration. 

A fast  Fourier  transform,  such  as  given  in  Ref.  [26],  is  a key  element  of  an 
efficient  implementation  of  transformations  using  Eqs.  (15a),  (15b),  or  (18). 
The  program  listing  contained  within  pages  202  to  240  of  Ref.  [26]  uses  fast 
Fourier  techniques,  and  makes  a good  starting  point. 

The  LINPACK  package  available  from  the  Society  for  Industrial  and  Applied 
Mathematics  (SIAM),  Ref.  [53],  contains  linear  equation  Bolvers  that  can  be 
useful  in  mapping  operations.  If  one  uses  a near-circle  to  circle  transform 
at  unevenly  spaced  points  in  the  circle  plane  (namely  at  points  corresponding 
to  the  input  points  in  the  physical  plane),  fast  Fourier  techniques  cannot  be 
used.  By  using  LINPACK  routines  on  a fast  computer  (i.e:  IBM  3081)  moderate 

computation  times  (about  one  minute  of  CPU  time)  are  required  for  a representa- 
tion involving  one  or  two  hundred  points.  A typical  FFT  mapping  time  for  the 
same  nuntoer  of  points  is  two  seconds. 

Fast  Poisson  solvers  are  available  for  filling-in  the  grid,  as  mentioned 
previously.  A rather  versatile  fast  Poisson  solver  is  reported  in  Ref.  [42]. 

Good  graphics  routines  for  plotting  intermediate  contours  and  grids  have, 
in  practice,  been  found  to  be  the  most  important  element  of  the  process  of 
inventing  new  mappings,  or  new  combinations  of  mappings.  By  using  extensive 
graphics,  a defective  mapping  is  soon  discovered  with  little  effort.  This 
allows  a wider  range  of  functions  to  be  tried  (guessed)  in  the  given  amount  of 
time  allowed  for  completion  of  a project.  The  graphics  package  must  be  general 
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enough  to  plot  both  contours  and  grids  and  should  be  accessible  by  a simple 
subroutine  call.  It  is  not  unusual  to  plot  fifty  different  intermediate  results 
during  the  development  of  a mapping  sequence  such  as  that  in  Ref.  [16].  Of 
course,  each  plotting  call  is  sequentially  commented  out,  but  not  deleted, 
during  the  development  process.  Such  an  extensive  use  of  graphics  greatly 
simplifies  debugging,  but  would  be  unbearably  tedious  without  a simple  calling 
sequence.  Since  graphics  tends  to  be  tailored  to  a particular  computer 
installation,  a sufficiently  versatile  and  easy  to  use  package  often  must  :e 
written  locally  and  thus  is  not  usually  available  to  the  public. 


CLOSING  REMARKS 

Modern  finite  volume  flow  solvers  do  not  require  an  orthogonal  grid,  but  a 
near-orthogonal  grid  is  usually  beneficial.  Often  a simple  conformal  trans- 
formation, followed  by  independent  algebraic  transformations,  can  be  used  to 
generate  such  a near-orthogonal  grid  with  little  difficulty  to  be  expected  in 
a computer  program  implementation.  This  is  a good  way  to  get  involved  with 
conformal  mapping. 

The  conformal  mapping  of  a contour  onto  a canonical  contour  is  far  easier 
to  accomplish  using  a one-step  technique  based  on  Eq.  (18),  as  opposed  to 
classically  mapping  a contour  onto  a near-circle  and  then  mapping  the  near- 
circle onto  a circle.  This  is  especially  true  since  Ref.  [26]  contains  a one- 
step  conformal  mapping  computer  program  which  is  easily  modified  to  map  new 
geometries.  This  one-step  mapping  technique  is  faster,  simpler,  more  resistant 
to  root  selection  problems,  and  more  stable  than  the  classical  technique. 

It  is  not  necessary  to  generate  the  grid  using  the  inverse  of  the  transform 
employed  to  map  the  contour  to  a canonical  contour.  In  fact,  it  would  seem 
that  use  of  a fast  Poisson  solver  to  generate  the  grid  using  the  known  boundary 
correspondence  offers  a flexibility,  simplicity,  and  economy  that  may  not  be 
surpassed  by  other  methods. 
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APPENDIX  A 

PRODUCT  MAPPING 

Consider  the  product  mapping 

f(Z)-f(Z)  « g(C)'g(C)  - w - <p  + iiji  . 

The  real  axis  in  the  Z plane  can  be  specified  by  the  relation  Z « Z,  so  on 
the  real  axis 

f (Z) • f (Z)  - f (Z) • f (Z) 

- f(z)’Kz) 

■ pure  real  ■ $ + ii|> 

Therefore,  the  real  axis  Z « Z maps  to  the  real  axis  iji  = 0 in  the  ui-plane. 
In  a similar  manner,  the  real  axis  4 “ C also  maps  to  the  real  axis  ij;  ■ 0 in 
the  a plane.  Thus  the  real  axis  Z ■ Z maps  to  the  real  axis  ; » 4. 

RATIO  MAPPING 

Consider  the  ratio  mapping 

im. « Aisi „ « „ 

f(Z)  I(Z) 


135 


APPENDIX  A (continued) 


The  real  axis  in  the  Z plane  is  specified  by  Z « Z,  so  on  the  real  axis, 


in 


f(Z) 

?(Z) 


in 

in 


f(Z) 

W) 

f(Z) 

rm 


■ pure  imaginary  ■ $ + it|i 


Therefore,  the  real  axis  Z « Z maps  to  the  imaginary  axis  ij>  » 0 in  the  u 
plane.  In  a similar  manner,  the  real  axis  ? « £ also  maps  to  the  Imaginary 
axis  4>  ■ 0 in  the  id  plane.  Thus  the  real  axis  Z « Z maps  to  the  real  axis 

t - l. 


